Charge dynamics of doped holes in high Tc cuprates - A clue from optical conductivity 
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The charge dynamics in weakly hole doped high temperature superconductors is studied in terms 
of the accurate numerical solution to a model of a single hole interacting with a quantum lattice 
in an antiferromagnetic background, and accurate far-infrared ellipsometry measurements. The 
experimentally observed two electronic bands in the infrared spectrum can be identified in terms 
of the interplay between the electron correlation and electron-phonon interaction resolving the long 
standing mystery of the mid-infrared band. 
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It is now widely recognized that the physics of dop- 
ing holes into a Mott insulator is the key concept to 
understand the high temperature superconductivity in 
cuprates An appealing scenario is that the spin sin- 
glet pairs already existing in the insulating antiferromag- 
net turn into the superconducting Cooper pairs when the 
doped holes introduce the charge degrees of freedom. On 
the other hand, it is also noted that the parent compound 
is an ionic insulator, where the polar electron-phonon in- 
teraction (EPI) plays an essential role, and it is expected 
that this strong EPI continues to be of quite vital im- 
portance even at finite doping. Therefore, the quantum 
dynamics of the doped holes is essentially influenced by 
both the magnetic fluctuations and quantum phonons. 

It has been recognized that the charge dynamics is de- 
termined not by the large Fermi surface but by the doped 
holes in the underdoped region [^. Therefore, it is a 
reasonable approach to consider the charge-current dy- 
namics of the holes (not electrons) interacting with the 
quantum phonons and magnons simultaneously to ana- 
lyze the infrared optical conductivity (OC). The basic 
features of the observed d.c. conductivity and OC fol- 
low Ref. 0, i, H. The OC cr(w) in undoped material 
reveals the charge transfer band at w = 1.5eV between 
the p-orbitals of oxygen and d-orbitals of copper. With 
doping, a low energy part of f7(w) develops revealing the 
dynamics of the doped holes as a function of the fre- 
quency uj. In particular, the Drude weight is shown to 
be proportional to x even in the absence of antiferro- 
magnetic long range order, while the relaxation rate 1/t 
is proportional to the temperature T. Therefore, this 
dependence of the Drude weight finds a natural expla- 
nation as directly reflecting the hole concentration, as 



mentioned above, while the T-dependence is due to the 
lifetime of the holes. The higher frequency cr(w), on the 
other hand, has not been well-understood. Especially 
the mid-infrared (MIR) peak with dependent on doping 
energy at around cjmir — 0.5eV is still controversial [5J, 
with interpretations involving an tj-dependence of l/r, 
transitions between the Zhang-Rice singlet state to the 
upper Hubbard band, and magnon sidebands. 

For this problem, angle-resolved-photoemission spec- 
troscopy (ARPES) offers an important clue. ARPES in 
undoped parent compounds measures the spectral func- 
tion of a single hole left behind when an electron is kicked 
out from the sample by the incident light 0, [3| . The cor- 
responding theoretical analysis has pointed out the role 
of the composite polaronic effect due to electron-magnon 
and electron-phonon (el-ph) couplings [1, Therefore, 
it is expected that these two interactions are indispens- 
able to understand the infrared optical spectra as well. 

We study in this paper theoretically the OC of the 
single-hole doped into a Mott insulator described by the 
t- J-Holstein model, and its numerical solution in terms of 
the Diagrammatic Monte Carlo (DMC) simulation. Pre- 
viously, the OC of the t-J-Holstein model has been cal- 
culated by exact diagonalization of small clusters [lo| . 
in the non-crossing approximation NCA for both mag- 
netic and lattice variables [ij], and for the case of infi- 
nite dimension [l^ . Compared with these approximate 
methods, our DMC simulation provides more accurate 
solution for the infinite system without approximation 
associated with the phonon sector [l^l and for proper lat- 
tice geometry and dimension. The only approximation is 
the NCA for magnons which is shown to be sufficiently 
good for the parameters considered below [3,111]. These 
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FIG. 1: Comparison of typical OCs of different models in 2D 
with experimental data of heavily underdoped cuprates: (a) 
Holstein model at A = 0.44; (b) t-J model at J = 0.3; (c) 
t-J-Ifolstein model for J = 0.3 and A = 0.39; (d) in-plane 
OC of 1.5% hole doped (Eui_^Ca^)Ba2Cu306 at T = lOK. 
The energy dependence of the theoretical data is presented in 
wavenumbers assuming t = 0.3eV (leV=8065.5cm-^). The 
absolute value of the theoretical ai is evaluated using the 
experimental hopping distance a — 3.86^4 and bulk hole con- 
centration rih = 1.72 X 10"^^ cm""^. 



results are compared "with accurate far-infrared ellipsom- 
etry measurement as "well as -with previously published 
data. The infrared ellipsometry measurements have been 
performed with a home-built ellipsometer attached to a 
Bruker Fast-Fourier spectrometer at the IR beamline of 
the ANKA synchrotron at FZ Karlsruhe, D at 70-700 
cm~^ and with conventional light sources at 500-7000 
cm-i [16,]. 

Figure [1] summarizes our main results, where four pan- 
els for the infrared parts of cr(cij) are compared. Fig- 
ures H^a), (b) and (c) show the calculated a{uj) for 
the Holstein model (el-ph coupling only), t-J model (el- 
magnon only), and the t-J- Holstein model (el-ph and el- 
magnon couplings) , respectively, while Fig. [Ijd) presents 
the experimental observation. Neither the t-J model 
[IJb) , nor the Holstein model [D^a) bear resemblance to 
the experiment [IJd), while their combination [IJc) at 
least qualitatively reproduces the salient experimental 
features. The clear signature of the experiment is that 
there are two prominent electronic components, i.e., the 
so-called MIR band at around o^mir — 4600cm^^ and the 
lower energy one at around uj = lOOOcm"^ that is located 
just above the infrared active phonon modes which show 
up as sharp peaks below 800cm^^. The lower energy 
peak roughly corresponds to that seen in Fig.[TJc) due to 
the phonon side-band. However, the MIR peak does not 
correspond to that of t-J model which occurs at around 
ojt-3 = 2 J = 2000cm~^. Instead of that, according with 
experiment, it is shifted to higher energies. 



As it can be seen from the results below, the coupling 
to two kinds of bosonic excitations results in two sepa- 
rate peaks because of the essentially different nature of 
the electron-magnon and el-ph couplings 17| and signif- 
icantly different energy scales of the magnetic and lat- 
tice excitations which are involved in the optical tran- 
sitions. The magnons with large characteristic energy 
~ 2J are weakly bound to the hole. To the contrary, 
the phonons are adiabatic and the EPI is considerable. 
As shown below, the lower energy peak is the phonon 
sideband with the threshold at the phonon energy [l^ 
and the higher energy peak is the magnon sideband of 
the lower peak. The reason for the apparent shift of the 
2 J peak to higher energies is most evident in the strong- 
coupling limit where the Franck- Condon picture for opti- 
cal processes is valid [l3| and the fluctuations of energies 
of different lattice sites, with the characteristic scale of 
Franck-Condon energy, can be considered as being frozen. 
Then, the energy cost of the transition of the hole from 
the ground state to excited states of the t-J model with 
frozen lattice is the sum of the energy of the emitted 
magnon and the Franck-Condon energy. Hence, the two 
peaks in the OC are the consequence of the importance 
of both el-ph and magnetic interactions. The same infor- 
mation is encoded in a different way in the single-particle 
spectral function observed in ARPES, where the low en- 
ergy quasiparticle peak of t-J model is affected by el-ph 
interaction. This low energy peak is separated into the 
broad Franck-Condon peak mimicking the dispersion of 
the t-J model while the zero-phonon line with very small 
weight has almost no dispersion 

In the standard spin-wave approximation for the t-J 
model 14, l3|, the dispersionless hole Eq = const (anni- 
hilation operator is /ik) propagates in the magnon (anni- 
hilation operator is ak) bath 
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with magnon dispersion cjk — 2 J^l — 7^, where 7k 
{coskx + cosfcj,)/2. The hole is scattered by magnons 
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with the standard scattering vertex Afk.q [4[. 

The OC of the t-J model has been calculated by var- 
ious methods in numerous papers giving mutually con- 



sistent results, e.g. [20|, |2]J. To confirm the validity of 
our method, we first reproduce the results for different 
J/t ratios [22]. We find the well known peak at around 
u!t-3 = 2 J, the origin of which has not been settled yet. 
It is also difficult to identify this peak with the observed 
MIR peak since the energy of the latter is about 2 times 
higher than 2J. Another source for skepticism is the 
opposite doping dependence of the Wt-j and wmir ener- 
gies 2l|. Therefore, we conclude that the t-J model can 
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not explain the observed OC even at very small dopings. 
Before introducing the EPI, let us first provide an inter- 
pretation of the Wt-j = 2 J peak. This peak comes from 
the hole excitations within the coherent band of the t-J 
model from the ground state at (7r/2,7r/2) to the neigh- 
borhood of (— 7r/2,7r/2) point, assisted by emission of 
single magnon with energy « 2 J and a momentum (tt, 0) 
(Fig. [2]). To prove this point we computed the OC by dis- 
entangling magnon and hole terms (exact at the lowest 
order in t) and then calculating the correlation function 
involving the hole operators in the subspace with one 
magnon (See upper panel in fig. [S]) 
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Here 



■0^^'' ) is the lowest eigenstate in the subspace with 



one magnon with energy 
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fco = (7r/2,7r/2) and 
23|. By direct inspection 
of the sum over q in Eq. [3] we found that the main con- 
tribution to the OC comes from the transfer of magnon 
momentums around (tt, 0). It can be traced out from the 
q-dependence of the vertex |C(ko — q, ko)p. 
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FIG. 2: (color online) Electronic transition (solid arrow) and 
emitted magnon (dashed arrow) responsible for the 2J peak 
in the OC of the t-J model. 

We now turn to the EPI added to the t-J model, the 
importance of which has been already established as men- 
tioned above. When the model is updated to the t- 
J-Holstein model, the hole interacts with dispcrsionless 
(frequency H. = const) optical phonons by short range 
coupling 7 
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In units of < = 1 we parametrize the dimensionless EPI 
constant as A = 7^/(4iri), choosing the value J/t — 0.3, 
and setting the phonon frequency O = 0.1 [il ]. 

Figure [3] shows the effect of EPI on the OC of a single 
hole in t-J-Holstein model. At weak EPI, an absorp- 
tion starts and shows a peak right above the phonon 
frequency. This apparent two-peak structure of the MIR 
response of underdoped cuprates can be tacitly discerned 




A,=0.441 
co/t 2 



FIG. 3: (color online) OC of a single hole in the t-J-Holstein 
model at J/t = 0.3 with various EPI coupling constants A. 
The vertical dashed line at uj/t — 0.1 indicates the phonon 
frequency. The dash dotted line in the upper panel is the 
result of Eq. ((3|. 



from many previous measurement (Fig. 3 in [25| and 
Fig. 9 in [2]) and is clearly seen from the low-temperature 
in-plane OC of 1.5% hole doped Eui_a;Caa;Ba2Cu306 
measured by ellipsometry (Fig. [Ijd)). This low energy 
EPI-mediated peak stays close to phonon energy up to 
the self-trapping transition point which is, for given pa- 
rameters of the model, is located at A « 0.4 Q. Indeed, 
according to the dependence of the dominant 2J con- 
tribution and the low energy peak on A (Fig. [3^), the 
transition from the weak- to the strong-coupling regime 
occurs at this coupling strength. We note that both the 
t-J model and polaron physics is crucial to explain the 
very existence of the two-peak structure of the OC. 

To understand the nature of the low energy peak in- 
duced by EPI, we did a calculation of OC for the Hol- 
stein model without hole-magnon interaction with re- 
duced transfer t = OAt mimicking the mass enhancement, 
which reproduces the self-trapping point of t-J-Holstein 
model with t — 1 As seen in Fig. 2)3, this effective Hol- 
stein model reproduces remarkably well the shape of the 
low energy feature of OC for the t-J-Holstein model. 

Since the effective EPI decreases with do ping [261 127| 
the reason of the experimentally observed [sl. I28l.l23| MIR 
mode softening is the change of the EPI. Comparing the 
position of the MIR mode with results of the t-J-Holstein 
model we can give a rough estimate of the renormaliza- 
tion of the effective EPI with doping. First, since the 
self-trapping point of the realistic extended tt't"-J model 
is A" * "-^ « 0.6 ^Sfli], we scale the EPI strength in the 
Fig. [4^ as A — > 1.5 A. Second, since the quadratic de- 
pendence of the energy scales of the OC response on A 
is a known property of the strong coupling regime [3] , 
we extrapolate the OC data to larger EPI couphngs. Us- 
ing a quadratic extrapolation of the experimental data of 
the MIR peak in YBCO [1, ^ and LSCO t29|] to zero 
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FIG. 4: (color online) (a) EPI coupling dependence of the 
energies of the dominating peak (solid line with circles) and 
EPI-mediated feature (dashed line with squares) of the OC. 
The energies of dominating phonon peak of the OC of the 
Holstein model at t = 1 (dashed line) and those at t — 0.4 
(solid line), (b) OC of the t-J-Holstein model (solid line) 
and OC of the effective Holstein model with t = 0.4 (dashed 
line) at the same coupling A — 0.1. Ratio of the effective 
EPI constant at doping x (or real in-plane concentration p 
for YBCO [sj]) to that at zero doping estimated from the 
MIR peak position for (c) LSCO and (d) YBCO mapped on 
the phase diagrams. 

dopings, we arrive at the result in Fig. \^-d which is 
in agreement with [26]. Moreover, since the analysis of 
the ARPES in undoped LSCO gives A - 1 the data 
in Fig. Ilt-d give absolute values of A'"^*^'-*(x). Finally, 
our result for ^lsco ^^j) jg quantitative agreement with 
the values obtained from the "kink" angle in ARPES on 
LSCO [i^. Figures IH^-d strongly suggest that supercon- 
ductivity appears after the effective EPI decreases from 
strong to weak coupling, which liberates the coherent mo- 
tion of the doped holes. The x-axis for YBCO should be 
translated to the hole doping concentration in the Cu02- 
planes p, which makes the two phase diagrams Fig. [D^c) 
and (d) look almost similar. Therefore, the behavior of 
the effective EPI and phase diagram seem almost univer- 
sal in high-Tc cuprates. However, the discussion above 
is restricted to the polaronic effect for the holes and not 
for the quasiparticles forming the large Fermi surface. 
The contribution of the EPI to the pairing hence is not 
excluded by the present analysis. 
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